Entanglement in Anderson Nanoclusters 
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We investigate the two-particle spin entanglement in magnetic nanoclusters described by the 
periodic Anderson model. An entanglement phase diagram is obtained, providing a novel perspective 
on a central property of magnetic nanoclusters, namely the temperature dependent competition 
between local Kondo screening and nonlocal Ruderman-Kittel-Kasuya-Yoshida spin ordering. We 
find that multiparticle entangled states are present for finite magnetic field as well as in the mixed 
valence regime and away from half filling. Our results emphasize the role of charge fluctuations. 

PACS numbers: 3.67.Mn, 71.10.Fd, 03.67.-a, 75.75.+a 
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In the last decade, solid state systems have emerged 
as a promising stage for quantum information process- 
ing, due to the prospect of scalability and integrability 
with conventional electronics. A successful realization 
of solid state quantum information processing requires a 
detailed control of the quantum mechanical properties of 
the system. In this respect, a key property is entangle- 
ment, or quantum mechanical correlations, between the 
individual quantum bits; entanglement plays the role of 
basic resource for a large number of quantum informa- 
tion schemes pj. In nanoscale solid state systems, inter- 
nal degrees of freedom of individual electrons, like spin or 
orbital states, are natural candidates for quantum bits. 
This offers considerable scope for studies of electronic en- 
tanglement in nanosystems. 

For spin entanglement, of particular interest are nan- 
oclusters with magnetic impurities, realized experimen- 
tally e.g. by magnetic atoms on a surface Q, on a nan- 
otube jij] or by coupled quantum dots 0. Such nan- 
oclusters display intriguing spin properties as the Kondo 
effect and antiferromagnetism, similar to what occurs in 
extended systems with dense magnetic impurities. For 
extended systems a central feature, described by the 
Doniach phase diagram ||, is the competition between 
formation of local Kondo spin singlets and a nonlocal 
Ruderman-Kittel-Kasuya-Yoshida (RKKY) spin order- 
ing. A model which captures such behavior is the pe- 
riodic Anderson model, PAM, a lattice of localized lev- 
els (with strong local repulsion) which hybridize with a 
conduction band. The PAM has been investigated in- 
tensively during the last decades |(| in connection with 
heavy fermion physics, non-fermi liquid behavior, etc 0. 

To date, no studies are however available for entan- 
lement in the PAM, only the related single [8j and two 
impurity Anderson models and the simplified Kondo 
necklace model ^3 have been considered. There are 
also investi gati ons of entanglement in various spin clus- 
ter models [ill H^ . None of these models do however 
correctly capture the interplay of spin correlations and 
charge dynamics, characteristic of the PAM 0. Only 
very recently the interest turned to systems where charge 
dynamics is important • These investigations how- 



ever only concerned the Hubbard model, focusing on en- 
tanglement in extended systems, at quantum phase tran- 
sitions [77) . An investigation of the entanglement in nan- 
oclusters described by the PAM is thus highly desirable. 

In this work we perform such an investigation. Us- 
ing exact diagonalization methods, ED, we study two- 
particle spin entanglement in the ground state and at 
finite temperatures. Both the Kondo- (well localized 
moments) and the mixed valence- (charge fluctuations) 
regimes are investigated. The effect of magnetic field and 
different electron fillings are considered. We find that: 
(i) The entanglement is governed by the hybridization, 
i.e. by the coupling between the conduction (c-) elec- 
trons and the local moments due to the f-electrons (see 
Fig. nj. At small hybridization the entanglement is pre- 
dominantly nonlocal, between f-electrons, while at large 
hybridization it is local, between f- and c-electrons. The 
local fc- entanglement survives at much higher tempera- 
tures than does the nonlocal ff- entanglement, (ii) These 
properties can be represented by an entanglement phase 
diagram, providing a view on the Kondo- RKKY competi- 
tion complementary to the Doniach phase diagram, (iii) 




FIG. 1: The concurrences Cg ] and C£ J as a function of 
k B T/t and V/t for L = 6, U = -2E { = 6t, B = and half 
filling. Projected at the bottom is the Doniach, DPD, and the 
corresponding entanglement phase diagram, EPD (see text). 
The L — 6 cluster is also shown. 
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In the mixed valence regime, at finite magnetic field and 
away from half filling we find nontrivial multiparticle en- 
tangled states, (iv) Charge dynamics plays an essential 
role in obtaining the results (i) to (iii). Our results illus- 
trate the richness of entanglement properties of the PAM 
and invite to further investigations of additional regions 
in parameter space, disorder and geometry effects, etc. 

Model and theory We study ring shaped clusters with 
L = 2 to 6 sites and two orbitals per site, denoted c and 
f (see Fig. \l]ior L — 6). Each cluster has N e electrons, 
with < N e < AL, The PAM cluster Hamiltonian is 

H = -* E (4<y* + h.c.) + uj2n(Ai + E fJ2 n L 

(i,j)a i i<y 

+ Vj2Ul^a + h.c.) +g^ B Bj2SL (1) 

ia io. 

with n\ a = fhfio and a =|, | , a = c, f. Here de- 
notes nearest neighbor, n.n., sites and the spin operator 
Sic = (Sf c ,Sf c ,Sf c ) = (1/2) £ CT , CT ' W ct^/ for c (simi- 
larly for f), with f = (t x , T y , t z ) a vector of Pauli matri- 
ces. The n.n. hopping between c-orbitals is t (t > 0), 
whilst V is the hybridization term between c and f- 
orbitals at the same site: both t and V can be taken 
real. The f-orbital onsite interaction strength is U, Ef 
is the energy of the f-orbitals and B is the magnitude of 
the uniform magnetic field. 

We focus on the reduced two particle spin density ma- 
trix p = p^a ) f° r orbitals a,/3 = f, c at sites i and j. 
Using the {ftt), ITI), lit), III)} basis, for e.g. a,(3 = f the 
density matrix elements are (f^f^, f ja „ f la ,„) where (..) 
denotes exact, many-body, thermal equilibrium averages. 
Due to translational invariance p depends only on \i— j\. 
Since [H, J2i a &ia\ = 0' eacn eigenstate has a well defined 
number of spin up and down electrons and only spin- 
conserving density matrix elements are nonzero. Con- 
sequently, p has in general five independent parameters: 



We point out that in contrast to the single site Fock 
space, or mode entanglement considered e.g. in 0, 113 . 
here the entanglement between two physical particles |2JJ| 
is considered. The concurrence C in Eq. is obtained 
from the reduced density matrix p which by definition 
determines any two-particle observable, as e.g. correla- 
tion functions. Thus, C(p) is a natural measure for the 
experimentally accessible two particle entanglement. 

Kondo-RKKY competition A central property of the 
PAM in the local moment (Kondo) regime is the com- 
petition between Kondo and RKKY correlations for the 
f-electrons. For macroscopic systems, such competition 
was described qualitatively by Doniach ||: for low tem- 
perature T and weak hybridization, V/t < 1, the local- 
ized f-electron spins are RKKY-ordered. Increasing V/t 
there is a cross over to local Kondo screening. Such com- 
petition occurs also in nanoclusters, where it is controlled 
both by V/t and the conduction level spacing [2l| • A de- 
scription in terms of a cluster Doniach diagram is possible 
|22j : A transition criterion was established by comparing 
spin correlators. Recently, similar results have been pre- 
sented for the 2D Kondo lattice model |23j ]. 

To investigate the entanglement in the Kondo-RKKY 
competition, we consider the symmetric case U = —2Ef 
(the ground state |^ g ) is a singlet |f|), with well localized 
f-electrons at half filling (N e = 2L), B = and T = 
0. Under these conditions the PAM can be mapped [24| 
onto the Kondo lattice model (KLM) characterized 
solely by J = 8V 2 /{Ut). In the KLM both the ff and fc 
density matrices are normalized and thus parametrized 
by a single parameter a. The concurrence is given by 
C = max{l— 6a, 0} and can be directly related to the spin 
correlator K = K^~ iY> = {S z la S z jl3 ) as C = max{-l/2 - 
6K, 0}, where K = a - 1/4. 

Preliminary insight in the Kondo-RKKY competition 
is gained from the KLM in the simple case L = 2, where 
l^/g) is obtained analytically. For the two limiting cases, 
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Note that typically tr(p) ^ 1, since the number of elec- 
trons at each orbital can vary from zero to two. When 
B = spin symmetry requires 16] the state to be charac- 
terized by only two free parameters, c and a, with a = a 
and b = b = a — c: the state is a Werner state |J3 • 

The entanglement is conveniently quantified via the 
concurrence C{p) = C^o - (p) = max{0, \Ai" - 
— \/^3 — V^ti} where the AiS are the real and positive 
eigenvalues, in decreasing order, of the matrix pp, where 
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(/iVi " /ii4t)( cT i T + C 2 T )( C H + c 2|)|0) J ■ 



(l/2)(/T T c T u - fl&Xfl^ ~ / 2 Vi T )l°> J -> °° 

(4) 

Thus, on increasing J from (small hybridization), the 
system evolves from a maximally entangled, singlet ff- 
state to a product state of cf-singlets (see Fig. At 
J w 1.22 we have C^ 1 = . Qualitatively this is 
similar to what occurs in the two-impurity Kondo |9| and 
Kondo necklace 10| models. 

For L > 3 (using numerical ED with U = — 2Ef = 6t), 
the physical picture for V < t is very different from the 
L = 2 case: The properties of \^> g ) depend strongly on 
size effects such as even-odd L parity and the density 
and configuration of the c-electron levels. This is clearly 

(i) 



( Ty ® T y )p*(T y ®T y ). This gives [3 for p in Eq. © illustrated by the drop of C£ from 1 for L = 2 to for 

L = 3, shown in Fig. [21 Interestingly, between L — 3 and 



C = 2 max{c- V^,0}, < C < 1 



(3) L — 6, increases monotonically with L from to 0.4. 
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FIG. 2: Left: Concurrence C ( a ' (for V = 0.2t) and C) c u) (for 
V — 1.5t) as a function of system size L. Inset: C$ } and c£ 0) 
as a function of J for the L = 2 KLM. Right: and C^ 2) 

and the total f-level occupation n/ as a function of £f/t in 
the mixed valence regime, V = O.lt, U = 6t and L = 4. The 
energy level scheme for V <C t is shown as inset: i? f (-Ef^) 
denotes the f-level in the mixed valence (Kondo) regime. 

This behavior follows from an increasing weight in |\I/g) of 
a superposition of a component with anti-ferromagnetic 
like order (not perfect for odd L) ~ | till ■••■) an d the 
same component with two neighboring spins flipped. 

We also note that, for L = 4 to 6, C s = while 
for L = 6 second next n.n. concurrence is again fi- 
nite, Cjp ~ 0.15. The nonlocal concurrences C^\ci 3 c 
with j > 1 are typically small or zero and will not be 
discussed further. It is important to note that |\& g ) al- 
ways contains some doubly occupied c-orbitals, making 
the result qualitatively different from the Kondo neck- 
lace model 0|- Also, an effective RKKY spin Hamil- 
tonian Q Jij^if ' &jf gives a very different \^> g ), 

and hence concurrence for most L. This provides strong 
evidence that the c-electron dynamics is essential for the 
obtained results. Both ff- and fc-entanglement decrease 
monotonically on increasing the temperature T, since un- 
entangled excited states get progressively populated. At 
ksT ~ O.lt (the typical energy scale of the low-lying 
spin excitations at V < t) ff- entanglement is suppressed, 
whilst fc-entanglement persists to much higher T. 

A global picture of the entanglement properties of the 
system is provided by an entanglement phase diagram, 
shown for L = 6 in Fig. ^ The phase boundary shows the 
crossover temperature for which Cs — Cf c . Importantly, 
for V < t the entanglement and the Doniach phase dia- 
gram (also in Fig. ^) have distinctly different crossover 
temperatures. Increasing T from the system thus 
passes three different phases: i) finite ff-entanglement 
and dominating ff-spin correlators | > l-Kfc | ; ii) 
zero entanglement and \Kg \ > | and iii) zero entan- 
glement and dominating fc-corr elations \K^ ] \ > \k£ } \. 
Moreover, there is a sharp cusp in the entanglement di- 
agram due to the nonanalyticity of the concurrence at 
C = 0. This illustrates the importance of the entangle- 
ment phase diagram |25| as a new tool for analyzing the 
Kondo-RKKY competition. 

Mixed valence regime. The effect of f-electron charge 
fluctuations on the entanglement manifests quite clearly 
when moving away from the Kondo regime. In the Kondo 



limit the f-level energies are well below the conduction 
band, |£f| ^> t, V and the occupation n/ = + is 
essentially one, due to the large value of U (= 6t) consid- 
ered. However, when Ef enters the conduction band, the 
regime of mixed valence, n/ starts to decrease. We stud- 
ied the mixed valence entanglement for —1 < Ef/t < 3, 
i.e. with nf < 1 (deep below the conduction band, 
Ej+U < —2t, the f-levels are doubly occupied and inert). 

We find that only for V < t there is additional entangl- 
ment away from the Kondo regime. In Fig. we plot C ff 
as a function of Ef/t for L = 4. When nf starts to drop, 
both n.n. and next n.n. entanglement increase, with a 
maximum at Ef w 0.2t and 1.7i, respectively. This can 
be understood from the manybody wavefunction around 
Ef ~ t, where the f-level occupation is roughly 1/2, i.e. 
on average two electrons occupy the f-levels. An analysis 
of the exact wavefunction shows that to a good approx- 
imation j^g) can be written as l^g) = \^f) <g> l^c), with 




a linear combination of f-electron singlets. The l^f) in 
Eq. JSJ gives a concurrence C^P = — 1/6, in rea- 
sonable agreement with the ED results in Fig. |2J The 
value 1/6 comes entirely from the low occupation of the 
f-levels, ti(p) = 1/6, since the singlet itself is maximally 
entangled. As in the Kondo regime, the entanglement is 
suppressed at fc^T ~ O.lt. The same result is found for 
all L, with the Ef/t interval of additional entanglement 
given by the density of states in the conduction band. 

Finite B-field Applying a magnetic field B, a naive ex- 
pectation would be that aligning more and more spins 
along the direction of the field would progressively sup- 
press the entanglement towards zero. It was however 
found in spin clusters that at strong magnetic fields, 
with all spins but one flipped, the state is ~ | J.TTT 
•■•) + I T1TT ■■■) + ■■, a strongly multiparticle entangled 
W-state |26|. This motivates an investigation the effect 
of B on the entanglement in the PAM. 

We plot (see Fig. 0) and as a function of 
B and V/t for L = 6. Since ^ la S z la ,H] = a finite 
B-field only shifts all the manybody states in energy an 
amount proportional to their total S z quantum number, 
eventually inducing a crossing of the B — levels and 
changing the ground state. This is clearly illustrated in 
Fig. |3| On increasing B, the total ground state spin 
S tot increments monotonically in steps of one, i.e. B in- 
duces successive flipping of the spins. While for V > t 
the fc-entanglement is suppressed on increasing B, the 
ff-entanglement survives for V < t up to 5 tot = 2. No- 
tably, for next and second next n.n. sites, entanglement 
actually increases with S tot and at S tot = 2 we have 
C { s j) w 0.3 for all j = 1,2,3 (Fig. 0). The analysis of 
the exact wavefunction shows that the ground state for 
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FIG. 3: Upper panels: The concurrence C l g ] (left) and c£ 0) 
(right) as a function of magnetic field B and hybridization 
V/t. Lower left: The ff-concurrence for j = 1,2,3 as a 
function of spin S tot . Lower right: The concurrence and 
as a function of hybridization V/t for N e = 8 electrons. 
In all panels, L = 6, E t = -U/2 = -3t and T — 0. 

V < t is given by |<F g ) = |%) <g> |<F C ), with 

i^> = E/iIMi°> (6) 

This is a W-state, just as was found in spin clusters. A 
magnetic field thus flips the f-spins before flipping any c- 
spin, a consequence of the much higher cost (t compared 
to O.Oli) in energy to flip a c-spin. fn the W-state all 
particles are mutually entangled with the same maximal, 
pairwise concurrence C = 1/3 |2*t| . in good agreement 
with the ED results in Fig. |3| The W-state is found 
only for L — 6; instead, for L = 3 to 5 there is negligible 
entanglement away from S tot — 0, an indication of a 
strong size dependence of this result. 

Away from half filling. For TV. ^ 2L, the PAM phase 
diagram is rather complicated [6|, in several parameter 
regimes the ground state is magnetic, S tot > 0. To in- 
vestigate the effect of the electron concentration on the 
entanglement, we considered fillings 2L — 1 > N e > L + 1 
for U = -2E f = 6t, B = and T = 0. We found that 
the overall behavior with ff- entanglement at V < t and 
fc- entanglement at V > t prevails away from half filling. 
For V > t, reducing N e monotonically suppresses the fc- 
entanglement, due to an incomplete local Kondo screen- 
ing. Interestingly, for V < t wc find finite ff-entanglement 
only for L = 4, 6 with N e = L + 2. Since « 1, this 
corresponds to two electrons in the conduction band. For 
this filling, the KLM ground state exhibits ff-spin di- and 
trimerization for L = 4 and 6 respectively |2Sj . We in- 
deed find that the many-body ground state is dominated 
by the term |'F) = j^f) ® l^c), with |*&f) a linear combi- 



nation of L dimerized terms, 

|*f) = 2- 1 [| 4ITT) + | TTU} + I ITT4) + 1 4TT4>] (7) 

for L = 4 and similar for L = 6. This state has the 
property that only sites opposite to each other in the ring 
show a finite entanglement, C^^ 2 ^ = 4/L, in reasonable 
agreement with the numerical results (see Fig. [3Jl- 

In conclusion, we have studied two-body entanglement 
in Anderson nanoclusters. We presented an entangle- 
ment phase diagram describing a generic entanglement 
scenario for systems with Kondo-RKKY competition. 
We also showed that Anderson nanoclusters exhibit mul- 
tiparticle entanglement depending on parameter regimes, 
electron filling or applied magnetic fields. More generally, 
our results give evidence that the interplay of charge and 
spin degrees of freedom must be taken into account to as- 
sess the entanglement behavior of nanoclusters with mag- 
netic impurities. We acknowledge useful discussions with 
C-O. Almbladh. This work was supported by the Swedish 
VR (P.S.) and EU 6th framework Network of Excellence 
NANOQUANTA (NMP4-CT-2004-500198) (C.V.). 
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